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Abstract 

We formulate correlation functions for a one-dimensional interacting spinless fermion 
model at finite temperature. By combination of a lattice path integral formulation for 
thermodynamics with the algebraic Bethe ansatz for fermion systems, the equal-time one- 
particle Green's function at arbitrary particle density is expressed as a multiple integral 
form. Our formula reproduces previously known results in the following three limits: the 
zero-temperature, the infinite-temperature and the free fermion limits. 

PACS numbers: 05.30.-d, 71.10.Fd, 02.30.1k 



1 Introduction 

The exact computation of correlation functions for strongly correlated quantum systems has 
been one of the major problems for years. Although, in general, this is exceedingly difficult to 
achieve, several analytical approaches especially in ID quantum integrable models have been 
provided to derive exact or manageable expressions of correlation functions. For instance, 
the low-energy behavior of correlation functions for gapless models can be systematically 
obtained by conformal field theory [1—4]. On the other hand, for systems with finite spectral 
gaps, the long-distance and -time asymptotics are investigated by (finite-temperature) form 
factor expansions (see [5,6] for recent developments). 

An alternative approach, which has been developed in these several years particularly 
for the spin- 1/2 XXZ chain, is to combine the algebraic Bethe ansatz [4] with solutions to 
the quantum inverse problem for local spin operators [7]. Using this, Kitanine et al derived 
multiple integral representations for zero-temperature correlation functions of the XXZ chain 
with an external field [8-10]. Their representations can be regarded as natural extensions 
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of the results based on the q- vertex operator approach [11-13], which is restricted to the 
zero-magnetic field case. One of the advantages of this method is that the formulation can be 
flexibly generalized to finite temperature and/or time dependent case [14-17] by combining a 
lattice path integral formulation. Furthermore, by considering a continuum limit of the XXZ 
chain, correlation functions of the ID boson system with delta function interaction can be 
obtained at finite temperature [18] (see also [19] for the zero-temperature case). 

Beyond spin systems, more recently we further extended the method to the calculation 
of correlation functions for fermion systems. By use of the fermionic i?-operator [20] acting 
directly on the fermionic Fock space, we have derived multiple integral representations of zero- 
temperature correlation functions for an interacting spinless fermion model with arbitrary 
particle density [21]. In this paper, we generalize the former results to the finite-temperature 
case by use of the quantum transfer matrix technique utilizing a concept of path integral 
[22]. Especially considered here is the equal-time one-particle Green's function. Our formula 
agrees with previously known results in the following three limits: the zero-temperature, 
infinite-temperature and the free fermion limits. 

The layout of the paper is as follows. In the next section, we review the quantum transfer 
matrix method for the spinless fermion model, and express the correlation function in terms 
of matrix elements of the monodromy operator. In section 3, we present the key ingredients 
of the computation for the correlation function. The multiple integral representation for the 
equal-time one-particle Green's function at finite temperature with arbitrary particle density 
is summarized in the main theorem. In section 4, the three special limits are evaluated. 
Section 5 is devoted to a brief discussion. The detailed derivation of the multiple integral 
form is deferred to the appendix. 



2 Spinless fermion model 

In this section, the thermodynamics of an interacting spinless fermion model is formulated 
by the quantum transfer matrix method [22]. The two-point correlation functions at finite 
temperature are expressed in terms of matrix elements of the monodromy operator. 



2.1 Fermionic i?-operator 

The Hamiltonian of the interacting spinless fermion model on a ID periodic lattice with L 
sites is defined as 

H ° = f E {sVi + 4+i c ^ + 2A (Q - n i) Q - - i) } > ( 2J ) 

where cj and Cj are the fermionic creation and annihilation operators at the jth site, respec- 
tively, satisfying the canonical anti-commutation relations. Here t and A are real constants 
characterizing the nature of the ground state, and fj, c denotes the chemical potential coupling 
to the density operator rij = c^Cj. 
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The underlying integrability of the model (2.1) can be seen by introducing the fermionic 
i?-operator defined as 



R,.(\) = l-n-- nj + gh(A ~ v) (nj + n 3 - 2n I% ) + -^—-^d Cj + cjcf), (2.2) 

which acts on Vj <g) s Vj . Here Vk is a two-dimensional fermion Fock space whose normalized 
orthogonal basis is given by |0)fc and jl)^ := cjjO)^, where Cfc|0)fc = and (g) s denotes the super 
tensor product. Identifying the above Fock space Vj (respectively Vj) with the quantum space 
TLj (respectively the auxiliary space TCj), we define the monodromy operator 7? (A) acting 
on the space Vj<g) s (Vi (g> s Vi (g> s . . . (g> s Vl) as 

Tf(\)=R lL (\)...R l2 (\)R ll (\). 

Since the fermionic /^-operator (2.2) satisfies the Yang-Baxter equation [20] 

-Rl2(Al — A2)i?13(Al)i?23(A2) = ^23(^2)^13(^1)^12(^1 — A2), 

the transfer operator 

T R (A) = Str T 7; R (A) = y<0|r R (A)|0> y - j(l\T R (X)\l)j 

constitutes a commuting family: [Tr(A), Tr(^u)] = 0, where the dual fermion Fock space is 
spanned by fc(0| and with := k(0\ck and fc(0|cj. = 0. The Hamiltonian Hq (2.1) is 
expressed in terms of the logarithmic derivative of the transfer operator Tr(A): 



H = tsh( V )-^ In T R (A) 



, A = chr/. 

A=o 



This relation yields 

Tr(A) = T R (0) (l + j^H + 0(A 2 )) . (2.3) 

For later use, let us define another type of transfer operator Tr(A) [22]: 
T R (A) = StviR lL (-X) . . .^ 2 (-A)%(-A)] 

with 

R- (A) = R st - j (A) = 1 - n 7 - rij + - f t hA J n T + n,- - 2n 7 n,-) (cU - c 7 c 7 ), (2.4) 

v y 1 ' 1 3 sh(X + r]) y 1 3 1 3 ' sh(A + r/) v ^ 3 lh v 7 

where stj denotes the supertranspose with respect to the jth space. Note that Tr(0) (Tr(0)) 
is the right-shift (left-shift) operator, namely Tr(0)xj = Xj+iTr(0) (Tn(0)xj = Xj_iTr(0)) 
where Xj = Cj,cJ, and hence T£ 1 (0) = Tr(0). Using this together with the expansion (2.3), 
one finds the statistical operator e~ H l T (T: temperature) is given by 



e -^/T = e Et=iMc(l-2n j )/(2T) lim 



:/ li( A):/„ |a-| 



A=0 

where the Trotter number N is assumed to be N £ 2N. Note here that we have set the 
Boltzmann constant to unity. 
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2.2 Correlation functions at finite temperature 

To derive multiple integral representations of finite temperature correlation functions, here 
we describe how the correlation functions can be expressed in terms of the transfer operator 
formalism. 

Let us consider a two-point correlation function at finite temperature T > 0: 

(O m+1 0\) = Tve -H/T ( m - *)> ( 2 - 6 ) 

where Oj is a local fermion operator. Inserting the formula (2.5) into the above, and using 
the fact that the i?-operator (2.2) or (2.4) is Grassmann even, one obtains 

_tv _ r Str T,...,FTri,...,L{T L (0) • • • (T m+1 (0)O m+1 )T m (0) • • • (T^oj)} 

( m+1 l} "i str r> ... >w ih > ...,L{r L (o)-r 1 (o)} 

where the operator Tj (A) acting in the space (Vj(g> s • • • <g) s V^) <8> s Vj is defined by 

T 3 {\) = e^-^y^R^-^R^ (a - A ) • • • %(-A)i? Tj (A - £ 

= A(A)(1 - nj) + B(A) Cj - + c]C(A) + D(X) nj . 
The Yang-Baxter equation and its modification 

R3l{— A2)-R32(— Ai)i?i2(Ai — A2) = -Rl2(Al — A 2 )i?3 2 ( — Ai)i?3i( — A 2 ) 

yield 

r 1 (A 2 )r 2 (A 1 ) J R 1 2(Ai - A 2 ) = Rn(Xi - A 2 )T 2 (Ai)T 1 (A 2 ), 
and therefore the quantum transfer matrix defined by 

T(X) = Tr j T j (X) = A(\) + D(X) (2.7) 
commutes for different spectral parameters: [T(A),T(//)] = 0. Thus (2.6) reduces to 

/<n ,. Str T ^{T i -™- 1 (0)Tr m+1 {T m+1 (0)O m+1 }T m " 1 (0)Tr 1 {T 1 (0)O 1 t }} 
(O m+1 Ol) = Jjm S^TriO) • 

Let us consider the thermodynamic limit L — > 00. Since the two limits L — ► 00 and 
N ^ co are interchangeable [23, 24], we can take the limit L — > 00 first. In addition, we find 
that the leading eigenvalue of the quantum transfer matrix T(0) (written as Ao(0)) is non- 
degenerate and separated from the next-leading eigenvalues by a finite gap even in the Trotter 
limit N — > 00. In the thermodynamic limit L — > 00, therefore, (2.6) can be written in terms of 
Ao(0) and the corresponding (normalized) eigenstate |^o) (note that Ao(A) := (^E r o|^(A)|\E f o))- 
Namely 

(0 n U = lim ^ol^ T n+i{T m+1 (O)0 m+1 }(A + J P)" t - 1 (O)'Ir 1 {r 1 (O)O;}|^o) 

^ m+1 11 N^oo A™ +1 (0) ' ^ ' ' 

In particular, for the equal-time one-particle Green's function (set Oj = Cj), which will mainly 
be considered in this paper, one obtains 

(c m+1 c\) = ( Cl cl +1 ) = - lim A - m - 1 (0)(* |C(0)(A + Dr-\0)B(0)\* o ). (2.9) 

N— >oo 

Note that we have used Trj{2}(0)cj} = -C(0) and Trj{^-(0)ct} = B(0). 
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2.3 Diagonalization of the quantum transfer matrix 

To evaluate the correlation function (2.8) (or (2.9)) actually, one must investigate the leading 
eigenvalues A o (0) and the corresponding eigenstates |^o) °f the quantum transfer matrix 
T(0). In this subsection, we present a general formula describing the eigenvalues and the 
eigenstates of T(A). The leading eigenvalue Ao(0) is expressed via the solution to a nonlinear 
integral equation. 

Let us define the reference state as 

\n) := |0> x ® |% ® ■ • • 10)^® |%. (2.10) 

s s s 

Obviously (2.10) is an eigenstate of T(A) (2.7): 

T(A) = (o(A) +d(A))|fi), A(\)\n) =a(X)\n), D(X)\Q) = d(\)\Q), 

where 



AT 



N 



shA 1 2 mc jv J sh(A - 



P\ 1 T 



Mc 
2T , 



a (\) = {— J —— } e #, d(A) = ("l)M » ) e 

In the framework of the algebraic Bethe ansatz, the vector |{A}) constructed by the multiple 
action of B(X) on \Q), namely |{A}) = nfli ^ s an eigenstate of T(A) if the complex 



parameters {Xj}^ =1 satisfy the Bethe ansatz equation: 



a ( A i) ( i \M TT sh ( A i ~ X k + r]) 



d{Xj) J-J sh(A 3 - - X k - 7]) ' 

The corresponding eigenvalue of T(A) is written as 



A(A> = «<A) n + n s^^. c, 12 ) 

The Bethe roots {A} characterizing the leading eigenvalue A o (0) are given by solutions to 
the Bethe ansatz equation (2.11) in the sector M = N/2. Then the Bethe ansatz equation 
(2.11) and the eigenvalue formula (2.12) are exactly the same as the spin-1/2 XXZ chain 
given by the Jordan- Wigner transformation, and hence we can directly utilize the method 
as in [25-27], which makes the analysis possible even in the Trotter limit N — > oo. Let us 
consider the following auxiliary function 

„(A) = (-!)* ^n||^i±|, (2.13) 

which associates the Bethe roots {Aj}^^ with zeros of 1 + a(A). To study the analytical 
properties of this function, we need to know the distribution of the Bethe roots describing 
the leading eigenvalue. It has been numerically verified for a wide range of Trotter numbers 
that the roots are distributed inside C. For instance, at \i c = 0, the Bethe roots for the 
critical (off-critical) regime are located on the real (imaginary) axis. Thus we can safely 
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assume the above features hold for any Trotter numbers, from which the analytical properties 
of the auxiliary function are determined. Consequently one sees a(A) satisfies the following 
nonlinear integral equation: 

lnfl (A) = -^ + ^ln Sh(A + ^ )sh(A -^-/^ ^)ln(l + aH) 

T 2 sh(A)sh(A- 4+7?) JC 2TTish(X-LJ + 7 1 )sh(X-LO- V ) 

Here the contour C is taken, for instance, as a rectangular contour whose edges are parallel 
to the real axis at ±7ri/2 (respectively ±r//2) and are parallel to the imaginary axis at ±77/2 
(respectively ±00) for the off-critical regime A = chr/ > 1 (respectively for the critical regime 
< A = chr/ < 1) (see figure 1 for a pictorial definition). In (2.14) the Trotter limit N — > 00 
can be taken analytically: 

ln (x) = _/fc _ tsh 2 (7?) _ r du; sh(2r ? )ln(l + aH) 

{> T Tsh(A)sh(A + ??) Jc 2vri sh(A - u + if) sh(A - u - rj) ' l ' ' 

For later convenience, we also introduce another auxiliary function o(A) = l/a(A) satisfying 
the following nonlinear integral equation in the limits N — > 00 [15]: 

in a(\) = *L - tsh2{v) 4- ( dU sh(27?)ln(l + aH) 

{ j T T sh(A) sh(A -v) Jc 2vri sh(A - w + 77) sh(A -u-ij)' { ' 

By the above auxiliary function o(A), the leading eigenvalue Ao(0) of the quantum transfer 
matrix T(0), which is related to the free energy density / by / = — TlnAo(O), is expressed 
as the following single integral form: 

1 a (o) = — 4- ( ^^MMi±^Ml _£i _ [ di0 sh ^ h ( 1 + °H) (9 171 
01 j 2T + Jc 2vri sh(a;) sh(a; + rj) 2T J c 2vri sh(w) sh(w - 7/) ' 1 J 

Differentiating (2.17) with respect to the chemical potential /j c , one obtains the particle 
density (rij) 

duj Tsh^a^oH _ /" dw Tsh(r/)a Mc o(Lj) 



/"du; Tsh^aH _ /" 
' ; J c 2 7 rish(a;)sh(a; + r ? )(l + a(w)) 7 C 



c 27ri sh(w) sh(w — 77) (1 + a(u;)) 



3 Multiple integral representation 

Along the line developed in [15], we can derive a multiple integral representing the equal-time 
one-particle Green's function for the spinless fermion model. Here we sketch briefly how to 
derive the multiple integral by presenting some crucial formulae to evaluate the action of 
the operator A + D and the resultant scalar product. These formulae are essentially the 
same with those for the zero-temperature case [21], since the commutation relations of the 
operators A, B, C and D are exactly the same with those for the zero-temperature case. 

First, it is convenient to introduce the following more general function $tv ({£}) instead 
of (2.9): 



(ttolC(ft) U?= 2 (A + fl)&)g(6n+i)|tto) 

irffAote) 



<M{£}) 

'i;= 

<{A}|C(&) YEUiA + D)(fc)B(Un)|{A}> 



(3.1) 



6 



(a) 



Im 



17T 

-~2- 



(b) 



Im 



Re 



[A 



-►Re 



Figure 1: The integration contours for the off-critical regime A > 1 (a) and for the critical 
regime < A < 1 (b). 

where {^j}™^ 1 is complex parameters located inside C. Note that {A} and |{A}) are, respec- 
tively, the Bethe roots and the eigenvector (not normalized), which characterize the leading 
eigenvalue Ao(0) (see the preceding section). The dual vector ({A}| is constructed by the 
multiple action of C(A) on the state ({\}\ = (^\Ylf=iC(Xj). It immediately follows 
that the one-particle Green's function (2.9) can be obtained by taking the homogeneous limit 
{£} -> and the Trotter limit N -> oo in (3.1): 

(ci4 +1 ) = hm lim<M{£})- (3.2) 

To evaluate the multiple action of the operator A + D on the state ({A}|C(£i), let us 
introduce the following proposition, which is originally proposed in the calculation of the 
correlation function for the spin- 1/2 XXZ chain [9]. 

Proposition 3.1 [21]. The action ofYl™ =1 (A + KD)(l;j) on a state {^\YljLiC(iJ,j) = ({^}\, 
for any sets of complex parameters {fij}^L 1 (not necessarily the Bethe roots), is written as 

m v 

(Mi + = E E Rn({t + }\{c}\{v + }\{v-}m + } u 

3=1 ™=0{ M }={^+}U{ M -} 

«}=« + }u{n 

where p = min(m, M), ({£ + } U = (0| YYj=i Ilfeii" ^(^k ) an< ^ coefficient R n 

is given by 

m—n i- n M—n 

Rn({C + }\{r}\{^}\{^})=Sn({C + }^ + }\{^}) II II ffa'tl) 

3=1 L k=l fc=l 

M—n 

(3.3) 

k=i fc=i 



7 



Here S n is defined as 



U] <k sh(^-^+) S h(e+-e fe + ) 



with 



M-n 



a=l 

M-n 



0=1 6=1 



sh(^+ - C 6 - ??) 



The functions /(A,//) and t(X, fi) appearing in (3.3) and (3.4) are, respectively, given by 

sh(A-^ + ??) shr/ 



sh(A — /x) ' sh(A — /x) sh(A — yu + 77) " 

Compared with that for the XXZ chain [9], some sign factors appear in the second term of 
(3.3) and (3.4), which originate from the fermionic nature of the present system. By setting 
k = 1 and applying the above formula to (3.1), one has 

$ m = _V V }|{A+}|{A-})({fi} U {A-}|g(Wi)|{A}> , . 

tt}=«+}u«-} 

|A+| = |?+|=n 

Here some new notations are adopted: 

(Al, . . . , Ajv/2+i) = (Ai, . . . , Ajv/2)£l)> (£l> • • ■ ,£m-l) = (£2, • • • ,£m)- 

Next we evaluate the action of £(£ m+ i) on U {A~}| by using the formula [21]: 

M M+l -f-fM ,/ . x 

<«i n^os^ + i) =(-i) M - 1 e ^.) n a^" M ? 

,■=1 j=i Ilfc=i sh{ni-n k ) 

M+i , , nlli 1 sh(^- - w + v) M+i 

x V - a[lll ' ] <«| TT Cfo), (3.6) 

where {Mjjjfji 1 are arbitrary complex numbers. One sees that the resulting equation consists 
of the ratio of scalar products such as ({£+} U {A~}|{A})/({A}|{A}), where {A} = {A+}™ =1 U 

{^J}f=i~ n and {£^}™ =1 € {^j}™^ 1 (see the appendix for detail). In fact, this quantity can 
be calculated by the following determinant representation of the scalar product. 
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Proposition 3.2 [21]. The scalar product between a Bethe state and an arbitrary state 



M M 

Sm(M|{a» = (n\ H c( H ) J] B(x 3 )\n) 

3=1 3=1 

can be expressed as follows: 

m(m-i) d(A,-)a(u,-) nfl- 1 sh(A 7 - — u k + rj) 

sm(mi{a» = 3 „m 3 : :\ / /; det M *(Mi{A}), 

where {Xj}jL l are Bethe roots, are arbitrary complex parameters. The MxM matrix 

*({//}|{A}) is defined by 

* i *(wi{A})=*(A J - >M *)-(-i) ^-^n^^- 

and detju denotes the determinant of an M x M matrix. 

Applying this, and using the same technique proposed in [15], one obtains the ratio of 
the scalar products ({£+} U {A~}|{A})/({A}|{A}): 



({g+}U{A-}|{A}) 
<{A}|{A}> 



n 

3=1 



q(#)(i + a(fr)) 
a(A+)a'(A+) 



>< n 

j,k=i 



nn 

j=l k=l 



7(a-,4 + ) " 

./(A-,A+) 
sh(A+ - A+ 



where the function G(A,£) satisfies the following linear integral equation 

dto sh(2?7) 



det n G(A+ £+), (3.7) 



G(\,Q=t($,\) + J 



l c 2-k'i sh(A — a; + ??) sh(A — u — rj) 1 + a(w) ' 

which can also be written in terms of o(A) as 

do; sh(2r/) 
c 27ri sh(A — u; + rj) sh(A — w — rj) 1 + a(o;) ' 



G(\,o = -t(\,0- J 



(3.8) 



(3.9) 



Applying all the steps described above, we find that (3.1) can be reduced to sums over 
the partitions of the sets {A} and {£}, and its summand consists of determinants of matrices 
constructed by functions of {A} and {£} (see (A. 3) for example). In fact, by using the 
technique as in [15], these sums can be transformed to multiple integrals on the canonical 
contour C, where the Trotter limit can be taken analytically. The derivation is straightforward 
but has a lot of steps, here we only write down the final result. Namely, the function $at({^}) 
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(3.1) is represented by the following multiple integral: 

m—l , 1 , m n+1 



(_i r /• sh(o-6-r/) 



^ n!(n + 1)! ~f[ 2 7ri b-(£j) sh((j - &n+i) 

dcjj sh(u;j - Cm+i) /" d^n+i /" do; ra+2 

C» jL^ 27ri (! + sh(w,- — ^ i — 77) J c 2vri(l + a(w„+i)) 7 C 2vri(l + a{u n+2 )) 



n 



n"i 1 [siiK+i-0 + '7)sii("»+2-0-'))] »',:(Ml{c}) 



n"=i [sh(cj„ + i - Uj + 77) sh(u;„ + 2 - Uj - f])} sh(a; n+ i - u n+2 + rj) 
x det„+iM- ({u;}|{C})det n+2 [G{uj, Ci), . . . , G(wj, Cn+i), £ m +i)] , (3.10) 

where 

A = l S h(^ - £j ± 7]) 

n"=i nLi sh (^i - wji, ± »/) n"=i n*=i s mg - c& ± »?) 

and M~({(j}|{C}) is an (n + 1) x {n + 1) matrix whose matrix elements are given by 

{+(, , AU +(f, , , \ Y\ n sh^q-cjj-T?) i-rn+1 sh(cjj-Ch-T?) f . ^ 

M _ = I t(Wj,Cfc) +HCfe,^)ll a =i Bh( Uj -u a -ri) 11&=1 sh(C 6 -^-r;) for j < n 



t($L,Cfc) for j= n+1 



o(A)(= l/a(A)) and G(A, £) satisfy the integral equations (2.14) and (3.8), respectively. C is 
the canonical contour defined as in figure 1, and T encircles £ and does not contain any other 
singularities. 

The one-particle Green's function can be obtained from the above expression by taking 
the limits {£} — > and N — > 00 (see (3.2)). The latter means to take o(A) as a function 
satisfying (2.15). We thus arrive at 

Theorem 3.1 The equal-time one-particle Green's function of the spinless fermion model at 
finite temperature has the following multiple integral representation, 

+ z. n!(n+1)! y rn+1 ii 27ri ^ sh(c .) ; 

f TT d^- / sh(o; i ) \ m /" du; ra+ i /" du n+2 

X i C n 1 = 1 27ri(l + a(uj)) \sh(uj - 77) y 7 C 2vri(l + oK +1 )) 7 C 2vri(l + a(u n+2 )) 

]l"=i [shK+i - Cj + rj) shK+2 ~ Q ~ *7)] W-({a;}|{C» 
n™=i [sh(w„+i - Uj + 77) sh(cj„ +2 - Uj - 77)] sh(w n+ i - u> n+2 + 7?) 

x det„ +1 M- ({a;}|{C» det„ +2 [G(ujXi), G(uj, Cn+i), G( Wj , 0)] , (3.12) 

where a(A) = l/(o(A)) and G(A,C) satisfy the integral equation (2.15) and (3.8), respectively. 
C is the canonical contour and T surrounds the point 0. 
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Using the identity 

1 K-^ = 1 - i — =TT> ( 3 - 13 ) 

l + a(w) l + a(cj) v ; 

we can convert the above multiple integral representation into another form. Namely, insert- 
ing the decomposition (3.13) into the part YVj=i 1/(1 + a ( u; j)) °f (3-10), and then performing 
the integrals over Q, we transform them to sums over the partition of the set {£}. Resumming 
the results in a similar way as in the appendix, we have 

n+1 



$ f muV ( ~ ir /" ff d ° sHCj ~ 6 + ^ 
" A. n!(n + 1)! 7 r „ +1 1 = 1 2vri b + (0) sh(C, - W) 



do^ b+(ujj) shjuj - Cm+i) /* dcj Tt+ i /• do; ra+2 



/" du; n+ i /" 
i c 2vri(l + a(cj n+ i)) 7 C 



X ^c™ n 27ri (l + fl(^i)) sh (^j - 6 + V) Jc 27ri (l + oK+i)) 7c 27ri (l + d(^n+2)) 

^ Il"=i [shK+i - Q + r?) sh(u; w+2 - Q - r,)] W+(M|{C}) 

n"=i [sh(w„ + i - + 77) sh(w„ +2 - Vj ~ v)\ sh(w n+ i - uJ n+2 + r/) 
x det„+iM+({a;}|{C})det n+ 2 [G{u jt &), . . . , G(wj, Cn+i), £ m+ i)] , 

where b+(w) and W+({u;}|{C}) are defined in (3.11). M + ({uj}\{(}) is an (n + 1) x (n + 1) 
matrix whose matrix elements are given by 

M+ = f *(a, + *k a) n^i SgpSS rei 1 fOT * * - 

3k \*(Cfc^i) for j = n+1 

Taking the homogeneous and the Trotter limits, we have another multiple integral represent- 
ing the one-particle Green's function. 

Corollary 3.1 The equal-time one-particle Green's function of the spinless fermion model 
at finite temperature has another multiple integral representation: 

"'- J (-if f Vj'dQ fsHQ + r,) 



(cic t } _y^ (-1)" f ffdo/^MO: 

X Jc" / = i 27ri (l + o(wj)) V sh K' + ") / 7c 27ri (l + o(w n+ i)) 7c 2vri(l + a(w n+2 )) 

E["=i [shK+i - + rj) shK+g - Q - ??)] W+(M|{Q) 
n™=i [ sll ( w n+i - w i + V) sh(w„ +2 - - rj)} sh(w n+ i - w„ +2 + n) 



x det n+ iM+(M|{C}) det„ +2 [G( Wj , Ci), • • • , G(uj, C„+i), 0)] , (3.14) 
J £i->o 

where a(A) = l/a(A) and G(A,C) satisfy the integral equations (2.16) and (3.9), respectively. 



4 Special cases 

In this section, we evaluate the three special cases of the multiple integral representation: the 
zero-temperature and the infinite-temperature and the free fermion limits. 
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4.1 Zero-temperature limit 

First let us consider the zero-temperature limit. Here we restrict ourselves on the off-critical 
case A > 0, and set 77 < as in [21]. The critical case < A < 1, of course, can be treated 
by just changing the definition of the integration contour C as in figure 1. 

Shifting the variables in (3.14) by Uj — > ujj — rj/2 and Q — > (j — rj/2, we deal with the 
integrals on the contour T v / 2 an d Co LiC^, where T v / 2 encircles the point rj/2; Co and C v /2 are 
defined as Co = [— 7ri/2, 7ri/2] and C v = [rj + 7ri/2,r? — 7ri/2], respectively. A close inspection 
of the auxiliary functions a(A) and o(A) for /j c > and 77 < at the zero-temperature limit 
T -» leads to 

1 t^o fl forAeZ 1 t^o fo for AG £ . . 

(4.1) 



l + o(A-2) forAG£ ' l + o(A-2) ^ for A G £ 

where £ = [—Q^ c , Q^ c ] and £ = (Co U C v ) \ £. Note that the Fermi point q^ c is an imaginary 
number (Img Mc > 0) depending on the chemical potential [x c . Substituting this into (3.9) 
and shifting the variables as above, one has 

H , 11 „ . f tW sh(2,)G(u; - 5 C-S) 



g(a-|,c-|) = -«(a,c) + /_ 



£ 2-7ri sh(A — u + rj) sh(A — u — rj) 



Comparing this with equation (2.24) in [21], one can identify G(\, Q with the density function 

G(A-5 C -5)=27rip(A,0. (4-2) 
Inserting both (4.2) and (4.1) into (3.14), we finally obtain 

m— 1 ^ „ rt+1 ]A . „ n+1 /■<//■ , ri 



£a<*o - e ^ i n £ /_ £ <^/ z «— n 

n=U 2 J = l J = l 

sh(^ - f)\ m nji 1 [siK+i - +>;)shK +2 - - 7,)] M'+({ W }|{c}) 



n 



t ^sh(uj + 2 ) J H™ =1 [sh(o; n+ i - u>j + rj) sh(a; n+2 - - rj)} sh(w n +i ~ + ??) 



x det n+ iM+({a;}[{C}) det n+2 \p(uj,Ci), ■ ■ ■ , p(e>j, Cn+i), p(wj, f?/ 2 )] • 

5 1— * 2 

The above representation completely agrees with equation (4.18) in [21]. 



4.2 Infinite-temperature limit 

Next we would like to deal with the infinite-temperature case T = 00, where the function 
3>jv({£}) (3.10) does not depend on the Trotter number N (note that o(A) = 1 and a(A) = 1). 
Therefore the integrals can be explicitly evaluated by just applying the residue theorem to 
the poles of integrand. The result reads 

lim (ci<4+i> = 0, 

1 —*oo 

as one expected. 
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4.3 Free fermion point 



Finally we investigate the representation (3.12) at the free fermion point A = 0. Set rj = 
7ri/2. Then the integral kernel in (2.15) and (3.8) becomes zero. Hence the two functions 
a(u>)(= l/o(ti;)) and G(oj,Q can be explicitly written as 



a(uj) = exp 



1 

T 



2it 



sh(2oj) 



G(u,Q = - 



sh(2(w-C))" 



Since Mjk = for j < n, one observes all the terms n > 1 vanish. Applying the decomposition 
1/(1 + a(u>i)) = 1 — 1/(1 + 0(0*1)), one finds that the integral including 1/(1 + cl(u>\)) is equal 
to zero since the integrand is antisymmetric with respect to uo\ and U2- After shifting the 
variables Q — > Cj + 7ri/4 and — > to, + vri/4, one obtains 



ch(o;i - C) ch(w 2 - C) 



sh(C 



1 



dcji 



do;o 



1 



sh(C + f i) J sh(2(C + f i)) J c > 2vri J c , 2vri 1 + a(uj 2 + f i 
1 1 



ch(wi-W2) Lsh(2(c^i -C))sh(2(w2 + fi)) sh(2(^ 2 - ()) sh(2(wi + § i)) 

where r_ 7ri / 4 surrounds the point ( = —iri/A; C = —Co U C_^ 2 ; Co = [—00,00]; = 
[— tti/2 — 00, — 7ri/2 + 00]. The integral with respect to u\ can be easily evaluated via the 
residue theorem applied to the poles at u> = — 7ri/4 and (. Then taking into account the pole 
outside the contour r_ 7ri / 4 i.e. at the point ( = u>2, we compute the integral with respect to 
(. It reads 



= -2(-iy 

+ 2(-l) r 



- f i) 
sh(w + I i) 



sh(2(^ + f i)) 



^ 2vri 1 + o(w + f i) 



dw 



1 



27ri 1 + ofu; 



sh(" - f i) 
sh(u; + f i) 

sh(w + f i 
~sh(w - |i 



1 



sh(2(o; + f i)) 

m ^ 

sh(2(w-fi)) - 



Changing the variable cosh(2u;) = 1/cosp (p £ [— vr/2, vr/2]) for the first term in the second 
equality, and cosh(2cj) = — 1/ cosp (p € [— 7r, —vr/2] U [vr/2, 7r]) for the second term, we obtain 



1 



( c i c m+i} = 



dpe 



imp 



1 



2vr 7_ 7r 1 + exp [-^ - f cosp] 2vr 7_ 7r 1 + exp + f cosp] 



dpe 



imp 



(4.3) 



We note that the above expression reproduces the well-known result (see [28] for example). 
Of course, (4.3) can also be derived by starting from (3.14). 



5 Discussion 

We have derived the multiple integral representation for the equal-time one-particle Green's 
function of the spinless fermion model at finite temperature. Unfortunately, the explicit 
evaluation of the multiple integrals still remains a difficult task, except for some special cases 
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considered here. Nevertheless, we believe that the method provided in this paper should be 
useful for the future study of the correlation functions of the fermionic systems. 

For instance, from (2.8) one sees that the long-distance behavior of the two-point corre- 
lation functions can be calculated by taking the ratio between the largest and the subleading 
eigenvalues of the quantum transfer matrix. For the one-particle Green's function, it reads 
(up to the sign) 

( c i c m+i) ~ 2A o cos(k F (m - 1)) exp 



with 



kp = Im 



In 



MO) 
A (0) 



Re 



In 



MO) 
Ao(0) 



<tfo|B(0)|*i) 



Ao(0) 

where Ai(0) is the leading eigenvalue for the sector M = N/2 — 1 (see (2.11) and (2.12)) and 
l^i) is the corresponding (normalized) eigenvector. In fact the finite temperature correlation 
length £ has already been calculated in [22]. The evaluation of the amplitude Aq by using 
(3.6) and Proposition 3.2 is quite important problem. 

It is also interesting to extend our result to the time dependent case. This is possible by 
combining the present method with the solution of the quantum inverse scattering problem 
for the operator Cj. It is evidently worth while to extract the long-distance and long-time 
behavior of the correlation functions at any finite temperatures. 



Acknowledgments 

This work is partially supported by Grants-in-Aid for Young Scientists (B) No. 17740248, 
Scientific Research (B) No. 18340112 and (C) No. 18540341 from the Ministry of Education, 
Culture, Sports, Science and Technology of Japan. 



Appendix. Derivation of multiple integral (3.10) 

We describe how the multiple integral representation (3.10) is derived. Applying the relation 
(3.6) to the term (Q\ \[ n j=1 C(§") Uk=i +1 ~ n C (K ) B (U+i) in the r.h.s of (3.5), we split 
<I>Ar({£}) into four parts according to whether the arguments of the functions a{x a ) and d{x d ) 
appearing in the resulting equation are Bethe roots {A} or inhomogeneous parameters {£}: 

$at({£}) =F 1 ({x a } e{x} \{x d } e{x} ) + F 2 {{x a } e{x} \{x d } em ) 

+ F 3 ({x a } e{(} \{x d } e{x} ) + F 4 ({x a } e{0 \{x d } e{0 ). (A.l) 
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First we consider the function F\ which can further be divided into two parts according to 
whether ft € {A+} or ft G {A"}: Fi = ^ ie{ A+ } + F £ l€ {A-}' where 

- 1 - ^({A-}!^})^^-}!^}) 



n=1 1A}={A+}U{A-} '=1 i'=l 



«}=tt + Mr} 

|A+|=|{+|=n 



({A}|{A})nr=i A o(0) 



2 ^ n 2 



xi?n(^ + }i{r}i{A + }i{A-})(oi n c(A-)n^;)c(^ + i)n B (^)i^) (a.2) 



with 



j=i j=i j=i 



d(K) Mv) n/=i /(a,-, a-) n -=i /(Af, tf) 

^ 1} ({A-}|{e + }) = _ - : , 

sh(A ; -Cm+ij 

a(A,7) sh^/^+x, A z 7) nf:r +1 /(Aj, A,7) n ;=i /(§"> A *~ 

^({A-}|U + }) = rr, jW ' : 

sh(^ m+ i - \, + 77) 



while F c ^ r T_i is 



F ^"l = ( " 1)T+1 S {5w £ {5 - } S I ({A}|{A})n™VA„te) 



«}=«+}u«-} 

|A+|=|{+|=n 



^n({^}i{r}i{A + }i{A-})^i n c(A-)n^;)c(ft, + i)ns(A J )in). 



J=l j=l j=l 



Here {A ± } denote {A ± } = {A*}^ 1 . 

Inserting the relations (2.12), (2.13), (3.3), (3.7) and a(A+) = -1 into (A.2) and shifting 
the variable n — > n + 1, we have 



JV _ AT _ 

m-2 2 n 2 Tt / i\ n n- 

p _ - \^ \^ \^ l -1 ) a l A z 

r 6e{A+i ~ 2^ Z> „i 



2 '"2 

6s(A+, = 5 {AW ^ u{A -> a '( A r)°'( A i~) n".i o'(a*) 

tt}=«+}u«-} 



|A+|=K+|-l=n 

^({A+}|{^})Z Tt ({A+}l{Q)y+({A+}|U+})X n ({A+}|{g+}) 
sh(A r - A7 +v)(l + a(ft)) nr=T~ 2 (l + 



(A.3) 
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where the functions y„({A+}|{£+}), F±({A+}| {£+}), X n ({A+}|{£+}) and Z n ({A+}|{£}) are 
defined as follows: 



with 



Y n ({\ + m + }) 



UU MA+) n"=i nS^(A+ - g - rj) sh(A+ - 4+ + V )} 

" n -i 1 n i ti - # + r?) n7,*=i sh (A+ - a+ - ^ 

x det n+1 M jk det n+2 G(Xj,^ k ) 



m— 1 



b ± (A) = n 



sh(A-^) 
i sh(A - £j ± ??) 



n 

J'=2 



sh(A-^) 
sMA-^ir?)' 



^+ £+m™ /( A ° ' A j + ) rr n+1 l^U^l 



/(A+Aj) ^"=i 

/(a£,|i) rrn+i /(a,g b + ) 



for j < n 



*(&.«.) + •<&>'<&. &) It, fgg IKS ^ fa i = n + 1 



and G(A, £) is the solution of the linear integral equation (3.8), where the variables {Aj}"^ 2 
and are, respectively, assigned as 



(Ai, . . . , X n +2) — (A^", • • • , A+, A ; , A ; , ), (£1, . . . ,^ n +2) — (Ci"> • • • >Ct+l' £rn+l)', 



n [i m lj mi sMe; - w n?=i bma; - a ± u) ' 
n{{ m }) "n^ 1 [ S h(A r -A;+r ?)S h(A--A;-,)]' 

in 

^«A + }|{e}) = n 



where |£ | = m — n — 2. Similarly, we obtain 



m— 1 



= 1 /(£,' ;A^) fc=1 f(^,Cj ) 



a(A" 



^'^^.hL-^S a'(A0a'(A,)n^a'(A;) 

tt}=tt+}u«-} 1+1 

|A+|=K+|=n 

. n t ({A + }i{r})^({A+}i{a)y+({A+}|{r})x„({A+}|{r}) 



sh(Ar _ A - +r?)n - T -i ( i +a(? - )) 



(A.4) 
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where (£, . . . ,£ +1 ) = , ■ ■ .,€+), = {^iJiT^, {£} = {O U {£"} and 

Y tsx+\\fr}) = UU b+(Aj+) UU n ^ [sh(A ^ + " ^ " V) ~ tk+ V)] 

m&) n "=i Mtf) n -ti - S + u) n ■ fc =i *ma+ - a+ - V ) 

X det„M_ 7 -fcdet n+ 2G ! (Aj,|fc), (|i, . . . ,£ n +2) = • • • ;£n+l>£m+l), 



Now we express the sums over partitions in (A. 3) and (A. 4) as multiple integrals. To 
this end, we would like to introduce the following useful formula. Let f(u>i, . . . , uo n ) be a 
function which is analytic on and inside the contour C, symmetric with respect to {^j}™ =1 , 
and zero when any two of its variables are the same. The poles of the function 1/(1 + a(uj)) 
inside C are simple poles at uj = Aj with residues l/a'(Aj), where {Aj}^ 2 are the Bethe roots 
characterizing the largest eigenvalue of the quantum transfer matrix. Hence one has 

|A+|=n 

The relation similar to the above also holds for a. 

First we apply (A. 5) to the partition for the set of the Bethe roots {A} in (A. 3). We 
see that the summand in (A. 3) has simple poles inside C at Xj = £7f. Since the inhomoge- 
neous parameters {£} can be chosen arbitrary values, we choose {£} such that the two sets 
of parameters {£} and {A} are distinguishable. Then there exists a simple closed contour 
surrounding the Bethe roots {A} but excluding {£}. Let C — T be such a contour, where T 
encircles {£}. Applying (A. 5) into (A. 3), one has 



m-2 



2 " 2 



F ^fx+} = Y Y Y Y ~\ I II 

K+|=n+l * 



2vri(l + a^)) 



(-l) ra a(Ar) 
a'(Aj-)c'(A,7) 



y n ({ w }i{e + })z n ({^}|{o)K + (V}iu + })A n ({^}|{e + }) 



( A - 6 ) 



sh(A- - A7 + r?)(l + a(6)) nr=T"ll + a«7)) 

By dividing the integrals, we transform the integrals along the contour C — T to those along 
the canonical contour C: 

n Tl -t Tl / \ n Tb k -J n k ~l 

UntHs-<)jLn£/.n%f". <-> 

where we have used the fact that the integrand in (A. 6) is symmetric with respect to {to}- 
Noting that, inside T, the integrand has simple poles at ujj = £7f , one can explicitly calculate 
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the integrals over T: 

k 



1% 



did, 



n—k+j 



2vri(l + a(iv n ^ k +j)) 



y„(V}|{e + })z„(v}|{e})y+(v}|{e + })x n ({^}|{e + }) 



{£++}u{£+-}={£+} 
l€+-|=* 



.7=1 1 + 11 



><n 



(5 



n—k 



m— n— 2 

n 



n-fc+l 



n 



a=l 6=1 /(£& >£ 



(A.8) 



By inserting (A.8) via (A. 7) into (A. 6), the integrals on the contour C — T can be transformed 
to those on the canonical contour C. 

The remaining task is the calculation of the sums over the partition of inhomogeneous 
parameters {£}. Resumming them by using the formula as in [15] 



Ixl \x~\ \x + \ 

E(- i ) fe e n V + "fixjMxj)] ui 1 - 

k=0 {x+}U{x-}={x}j=l 3=1 

\x+\=k 



X 3 



3=1 



\x\ 



and further expressing the sum over A ; (respectively A ; , ) as the integral over uj n+ i (respec- 
tively uJn+2) by (A. 5), one has 



rn-2 



{-iy 



-E E 

n=0 tt}=tt + }U«-} 3=1 
K+|=n+l 



7 c „ y [27ri(l + a(wj))J ic-r 



du;. 



n+l 



27ri(l + a(a; n +i)) 

dc n+2 K + ({^}|{e + })^-(v}|{g+})A w (v}iu + }) 

' c _ r 2vri(l + a(^ n+2 )) sh ( Wn+1 - u n+2 + r?) (£+)(! + 0(6 )) 

xdet„ +1 M«({u,} I {£+})det„ +2 [G(u;,-,£+ ),..., G(^-, G(^-, ^+i)], (A.9) 

where and the (n + 1) x (n + 1) matrix Afj^({w}|{^ + }) are, respectively, 

defined by 



UU mtl M", - ± v) MZt -uj±v) 
Ulk=i M"3 - w fc ± n) n • ,ti - 4 + ± v) 



M&\{»}\{ey) 



= < 



Hi i • 4- f C?+ r > TT n sh^-o^-??) i-rn+1 s\i{uj-£+-rj) 

^ 1 + )+(<:, ?+m n sh K-fr+ ? ?) rr n+1 sh (^"~^~^ 



for j < n 
for j = n + 1 
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The integrand of (A. 9) is a symmetric function with respect to {£ + } and vanishes when any 
two of them are the same. Thanks to this together with the fact that \/b{oj) has simple poles 
at uj = we can directly apply (A. 5) to (A. 9). Thus we arrive at 



m— 2 



i-iy 



■tie{\+} E„!( n+ i) 

n=0 v ; 



/ 

Jc-r 



dia 



'n+l 



/ 

Jc-r 



dC 



2vrib. 



(C J )Jy Cn f = i 1 [27ri(l + a(^))J 



du n+2 K+(M|{C})w„-(MI{C})^»(MI{C}) 



2vri(l + a(uj n+1 )) J c _ r 2vri(l + a(uj n+2 )) sh(cj n+ i - uj h+2 + v)( 1 + &(£i)) 
x det n+ iMW({w}|{C})det n+ 2[G(w j , Ci), • • • , G(uj, C„+i), G( Wj , U+i)], (A.10) 

where T = T — F^; surrounds the point £i but excludes {£}. 

Almost the same method is applied to F^^-i (A. 4) by considering the integrals over 

the contour C = C — T^ instead of C. Utilizing the transformation (A. 5), and resumming the 
resulting equation as in the case of F^ e ^ + y, one may have 



m— 1 



iie{\-} Yl Yl 



(-1) 



m— n— 1 



du>jb-(uj) 



'c» l>i(l + 



do; 



'n+l 



c _ r 2vri(l + a(cj n+ i)) 



I 



n=0 «}={J + }u«-} 

K+I=n 

d^ n+2 t/ n (v}i{r})K + (^}i{r})^({^}i{r})A„(v}i{r}) 



c _r 2vri(l + a(w n+2 )) 



sh(a; n+ i - cj„ + 2 + r/)b_(Ci) ]TLi b-(£ 7 : 



xdetjl4({u;}|{£ + })det^^^^ (A.ll) 

where the function U n ({uj}\{^ + }) and n x n matrix Mj/,({uj}\{^ + }) are, respectively, written 
as 



rr /r, ,Ti rT+i^ _ TT s H J ZSLZ '11 TT ail ^ fc ~ ^ 



+ »/)' 



sh(cj a - Wj - rj) ^y 1 sh(a; j - £ - rj) 



■+ 



a=l ^ - 



n 



Applying again the formula (A. 7) to the integration over C = C — T^, and noting that the 
sum over k in (A. 7) is restricted to k = and k = 1, we divide F^ e ^_y in (A.ll) into 

the following two parts: F c ^ { \-\ = r ~ , + F^ r ~ , > where ~ , is given by simply 
_ £i e l A ) &e{A-} 6e{A-}' 6e{A-} b J ^ J 

changing the contour C — > C in (A.ll), while ~ , is written as 

v ; ae{A-} 



m— 2 / 1 Xm „ n+l r i j. -i „ n 

F (i) = v / TT dCj / TT 

66{ A- } A, n!(n + i)l J fn+1 11 [2vrib_(0) J Jen f = \ 



dujjb-iujj) 



L 



duo 



n+l 



t 2mb-(C j )\J c »fJ i l2m(l + a(u j ))\ 

dco n+2 V+{{^}\{C})W-{{u}\{Q)X n {{u}\{Q) 



Ic-r 2vri(l + a(uj n+1 )) J C _ T 2vri(l + a(uj n+2 )) sh(u; n+ i - u n+2 + r?)(l + a(£i)) 
x det n+ iM^ ) ({a;}|{C})det n+2 [G(a; j , Ci), • • • , G(uj, ( n+1 ), G{uj, U+i)]- (A.12) 
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Note here that^we have shifted the variable n — > n + 1 and converted the sum over the 
partition for {£} into the integrals over F. The (n + 1) x (n + 1) matrix M-j^ ({u} \{(}) is 
defined as 



M^ } (M|{C}) 



IKVj^k) ^ z K^k,Uj) ll a= i sh^-ua-rj) llb=l B h{C b -Uj-ri) 



for j < n 



+(r, \ 4- hp, r, \ rr n M^a-ti+v) n n+1 sh (C6-6-^;) f nr _ „ , i 
+ Wllo=i shK-a-j?) llb=i sh(C-6+r;) J ~~ + 



In the next step, we would like to consider the sum i 7 ! = i* 1 , f= /T+\ + -f 1 < 



6e{A+}^ x 5ie{A-} 



and combine 



the three multiple integrals into one. First we deal with the sum F r ^ r r +1 + ~ . From 
(A. 10) and (A. 12), it immediately follows that 



dt^n+2 



f +pW -V ( ~ ir /" fff dC ^ 1 /" TT 

x f duj n+i r 

J c _ r 2m(l + a(u n +i)) J C - 



2vri(l + a(uj)) 



V-({co}\{C})W-({u;}\{C})X n ({co}\{C}) 



c _ r 27ri(l + a(uj n+ i)) J C _ T 2n\(l + a(u n+2 )) sh(w n+ i - w n+2 + rj) 

x det n+ iM^({w}|{C})det n+2 [G(wj, Ci), • • • , G(wj, Cn+i), G(wj,f m +i)], 

where the elements of the (n + 1) x (n + 1) matrix M~({u;}|{C}) are given by 



M _ = h(^,c*)+t(a,^)nLiS 



sh(i^ a -qjj-?7) y-jn+l sh(uj~Ci>-") 
'J " 



nn+1 siiy^ 3 -^b—<l) r ■ < 
_ _ , , 6=1 sh(C 6 -^.-r,) tor j < n 

for j = n + 1 



Changing the contour r — > T and combining it with F^ ~ , we obtain 

6e{A } 



m— 1 



Fi 



y ( " 1)T / n 



n+l 



do 



In 



L 



dui n 



+i 



L 



27rib_(0)J Jc^fJi L27ri(l + o(wj))J 

du n+2 y w -(M|{C})w-(M|{c})x w (M|{c}) 

/ c _ r 27ri(l + a(w n+ i)) J c _ r 27ri(l + a(w n+2 )) sh(w n+ i - w n+2 + r/) 

x det n+ iM- ({a;}|{C})det n+2 [G(a; j , Ci), . . . , G(wj, Cn+i), G^-.Cm+i)]- (A.13) 

The remaining contribution F 2 + F% + F4 in (A.l) can be absorbed into (A.13) by changing 
the integration contours for uJ n +i and uJ n +2 as C — T — > C. We thus finally arrive at 



m— 1 



(-1)' 



n=0 
dWn+1 



n\(n + 1)! Jpn+i 



n+l 

n 

j=i 

dw„+ 2 



do 



2vrib_(Ci 



n 

j=i 



do;jb_(o;j) 
2vri(l + a{ujjj) 



V-{{u}\{Q})W-{{u}\{Q)X n {{u}\{Q}) 



Ic 27ri(l + a(w n+ i)) 7 C 27ri(l + a(w n+2 )) sh(u; n +i - uj n+2 + r/) 

x detn+iM^ ({w}|{C})det n+ 2[G(wj, &), • • • , G(wj, Cn+i), G(w j , Cm+i)]- 
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